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I. INTRODUCTION 

The top quark is the heaviest particle so far discovered, with a mass close to the electroweak 
symmetry breaking scale, and closely related to various extensions of the standard model (SM). 
Thus, it provides an effective probe for the electroweak symmetry breaking mechanism and a test 
for the predictions of the SM through its production or decay. 

The production of the single top provides a good opportunity to study the charged weak current 
interactions of the top quark, e.g., the structure of the the Wtb vertex |fl|. Besides, it is an important 
background in many new physics searches at hadron colliders. However, due to the difficulties to 
discriminate its signature from the large background, it takes a long time after the discovery of 
top quark for the DO [Q] and CDF [|| collaborations at the Tevatron to observe the single top 
production. Recently, the ATLAS and CMS collaborations at the LHC have also measured the 
cross section of single top production at low integrated luminosities ||, ||] . 

Among the three production modes at hadron colliders, the t-channel is specially important 
because of its largest cross section at both the Tevatron and the LHC. This process has been ex- 
tensively studied, including the next-to-leading (NLO) QCD corrections based on the 2—7-2 leading 
order (LO) process, called the five-flavor (5F) scheme |^-13|. It is shown that the NLO corrections 



increase the LO cross section by about 9% and 5% at the Tevatron and LHC, respectively. In Ref. 
[14], the NLO calculation of the t-channel production based on the 2 — > 3 LO process, called the 



four-flavor (4F) scheme, was presented, which shows that the inclusive cross section in 4F scheme is 
smaller than in 5F scheme while the uncertainty in 4F scheme is larger than in 5F scheme. This is 
due to the fact that in the 5F scheme, the large logarithm of the form log(Q 2 /m|) due to the initial 
bottom quark is resummed into the bottom quark parton distribution functions (PDFs) and thus 
the scale dependence is significantly reduced. Besides, in the 5F scheme, the parton shower Monte 
Carlo simulation for the t-channel single top production was studied ]l5|-[r^|, and the threshold 



resummation for this process is carried out with the conventional resummation method [18-2C], 
where the partial next-to-next-to-next-to-order results are obtained by expanding the resummed 
cross sections to avoid the infrared singularities and ambiguities from prescription dependence. 

In this work, we investigate the resummation of the t-channel single top production in the 5F 
scheme using soft-collinear effective theory (SCET) [^-25]. SCET is developed to describe the 



behavior of the QCD interactions in collinear and soft regions with the short distance information 
encoded in the Wilson coefficients. It is very suitable to deal with the scattering processes with 
multiple scales. In the past ten years, SCET has proved very useful in high energy hard scatter- 
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ing processes. In general, these processes can be divided into two kinds, i.e., the time- like and 
space-like. The time-like processes produce a time-like particle in the intermediate or final state, 
including Drell-Yan production p^-p9[, Higgs boson production [27, |30|-34|, e + e~ annihilation to 



hadrons 1 35 -|39j[ , color-octet scalar production [40], direct top quark production via FCNC cou 



pling [41 1, and s-channel single top production j42|. The space-like processes involve a space-like 



particle in the intermediate state, such as deep-inelastic scattering [25, 43-45], direct photon pro- 



duction pf| and W (Z) boson production at large transverse momentum px p7[- Note that some 



processes are mixing of these two kinds, e.g., the top quark pair production |48-p0|. 

The threshold region can be easily defined for the time- like processes. It is usually defined 
as the limit z = m 2 /s — > 1, where m is the invariant mass of the time-like particle and s is the 
square of the center-of-mass energy. For the space-like processes, the threshold region is a little 
more subtle. The threshold region for the deep-inelastic scattering process is given by the Bjorken 
scaling variable x — > 1. For the direct photon production and W (Z) boson production at large 
Pt, the threshold region is approached when £4 = M 2 — > 0, where M is the mass of everything 



in the final state except the photon (W or Z). The t-channel single top production is a space-like 

x 



process involving four colored external particles. We define the threshold region as S4 = P\ — > 0, 



similar to the case of W or Z production at large pt, where P\ represents the mass squared of 
everything in the final state except the top quark. In this threshold region, the cross section can 
be factorized as 

a = H®J®S®f Pa ®fp b , (1) 

where H, J, S, fp are the hard function, jet function, soft function and PDF, respectively. The 
hard function incorporates the short distance contributions arising from virtual corrections. The 
jet function describes all collinear interactions inside the outgoing jet. The soft gluon effects coming 
from all colored particles are contained in the soft function. The PDF denotes the probability of 
finding a particular parton in the proton. 

The final states of the t-channel single top production at hadron colliders consists of a single 
top quark and a jet at the LO. Additional soft gluons can be emitted from the colored initial- 
and final-state particles, and collinear gluons can be emitted in the jet. These contributions are 
of higher orders in a s , but can be large in the threshold region. Besides, the hard part of this 
process receives a large correction since the usually chosen renormalization scale and the typical 
transferred momentum are rat and Q = \f—i, respectively, which yield a large logarithm of the 
form ln(m 2 /Q 2 ). Therefore, it is necessary to resum all these large logarithm to all orders. In the 
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SCET approach, the different scales in a process can be separated because the soft and collinear 



degrees are decoupled by the redefinition of the fields [24]. At each scale, one only need to deal with 
the relevant degrees of freedom. In this way, reliable perturbative expansions can be achieved easily, 
and the dependencies of the final results on the scales are well controlled by the renormalization 
group (RG) equations. As a result, the singular terms in the hard, jet and soft functions can be 
resummed conveniently. Furthermore, in numerical calculations, we find that for top quark px > 
50 GeV, the singular terms approximate the fixed-order calculations well, but for pr < 50 GeV, the 
singular terms do not dominate over the NLO corrections. This is understood because in the larger 
Pt region, the phase space for the additional emitted gluon is more constrained so that the main 
contribution comes from the soft gluon effects. Thus, we need to know an improved resummation 
prediction on the top quark transverse momentum distribution in the region of large pt, instead of 
the total cross section. Such t-channel top quark transverse momentum distribution is actually an 
observable which can be compared directly with the experimental results 1 , and is an important 
background in searching for new physics. For example, if there is an extra gauge boson W with a 
mass around 1 TeV and the standard model-like couplings, it is better to search this gauge boson 
in the tb final states than the two light jets final states because of the large dijet background from 
QCD processes. Moreover, one should concentrate on the events with large top quark px since the 
top quark from the decay of W usually has a large momentum. In this case, a precise knowledge 
of the t-channel top quark transverse momentum distribution at large pt in the SM is necessary. 

This paper is organized as follows. In Sec. O, we give a brief introduction to SCET. In Sec. ID, 
we analysis the kinematics of the t-channel single top production process and give the definition 



of the threshold region. In Sec. IV, we present the factorization and resummation formalism for 
the t-channel single top production in momentum space. In Sec. |V|, we calculate the hard and 
soft functions at NLO. Then, we study the scale independence of the final result analytically. In 



Sec. VI, we discuss the numerical results for t-channel top quark transverse momentum distribution 



at the Tevatron and the LHC. We conclude in Sec. VII. 



1 We have discussed with the ATLAS and DO experimentalists by email. They plan to provide such differential 
distributions in an update of the cross section measurement. 
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II. BRIEF INTRODUCTION TO SCET 

To describe collinear fields in SCET, it is convenient to define a lightlike vector n a = (1, n), n 2 = 
1. Any four-vector can be light-cone decomposed with respect to n M and n a = (1, — n) as 

with l + = n-l and l~ = n-l. The momentum of a collinear particle moving along the n^ direction 
has the following scaling 

p» = (p+,p-,p n± )~(\ 2 ,l,\), (3) 

while for a soft particle, the momentum scales as 

q ~ (A 2 , A 2 , A 2 ), (4) 

where A <C 1 is a small expansion parameter in SCET. E.g., for an energetic jet with invariant mass 
mj and energy Ej, A = mj/Ej. From the momentum scaling, one can see that the interaction 
between collinear fields of different directions rij and rij with rii • rij 3> A 2 will inevitably change 
the momentum scaling, thus is forbidden in SCET, but can be included as an external current in 
our computation. The soft fields, on the other hand, can interact with any collinear field without 
changing the scaling. 

In SCET, the n-collinear quark ip n and gluon field An can be written as 

Xn{x) = W n {x)£ n {x), £ n {x) = ^^„(x), 

KM = \w ] n iD» nl W n {x)\ , (5) 

where 

*Kl = K± + 9sKl (6) 

is the collinear covariant derivative and the label operator V is defined to project out the large 
momentum component of the collinear field, e.g., Vn£,n = P^Cn- Here we have split p into a sum of 
large label momentum and small residue momentum, 

p^=p^ + k /1 , With P^=P~Y +P n±- ( 7 ) 

The n-collinear Wilson line, 

/ r° \ 

(8) 



W n (x) = Pexp (^9s J dsn-A^(x + sn)t a ^j , 



6 



which describes the emission of arbitrary ra-collinear gluons from an n-collinear quark or gluon, is 
constructed to make the collinear fields as defined in Eq. (||) invariant under the collinear gauge 
transformation. The operator P is the path-ordered operator acting on the color generator t a . 
At the LO in A, only the n-A s component of soft gluons can interact with the n-collinear field. 



Such interaction is Eikonal and can be removed by a field redefinition [24|: 



where 



Xn(x) -> Y n (x-)xn(x), 

A» ± (x) -> y„(x_)^ ± (x)y n t(x_), (9) 

Y n (x) = P exp (ig s J ds n- A a s (x + sn)t°^ (10) 



for an incoming Wilson line [24, 51]. And for an out going Wilson line, it is defined as 

Y n {x) = Pexp (-ig s J dsn-A a s (x + sn)t a ^j. (11) 

The soft gluon fields are multipole-expanded around x_ to maintain a consistent power counting in 
A. For the interaction between the soft gluon fields and massive quark fields, there exist a similar 
timelike Wilson line [ ja2]j , for example, 

Y v (x) = P exp (ig s J dsv-A a s (x + sv)t a ^J . (12) 

After the field redefinition, the LO SCET lagrangian is factorized into a sum of different collinear 
sectors and soft sector, which do not interact with each other, 

£ S CET = ^4° l ) +£» + ••• • (13) 

The decoupling of soft gluons from collinear fields is crucial for deriving the factorization formula. 



III. ANALYSIS OF KINEMATICS 



In this section, we introduce the relevant kinematical variables needed in our analysis. As an 
example, we consider the subprocess 

<Pa) + b(p a ) -> t(q) + X, (14) 

whose Feynman diagram is shown in Fig. [l]. First of all, we define two lightlike vectors along the 
beam directions, n a and n^, which are related by n a = n^. Then we introduce initial collinear fields 
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ir 



FIG. 1: The Feynman diagram for the subprocess ub — » dt. 

along n a and n b to describe the collinear particles in the beam directions. In the center-of-mass 
frame of the hadronic collision, the momentum of the incoming hadrons can be written as 

PH = Ecm^, K = E CM f. (15) 

Here Eqm is the center-of-mass energy of the collider and we have neglected the masses of the 
hadrons. The momentum of the incoming partons, with a light-cone momentum fraction of the 
hadronic momentum, are 

p a = X a E C M~^, Pb = X b E C M-^-- (16) 

At the hadronic and partonic levels, the momentum conservation means 

P a + P b = q + Px, (17) 

and 

Pa+Pb = q + Px, (18) 

respectively, where q is the momenta of the top quark. We define the partonic jet with momentum 
px, which represents the set of all final state partons except the top quark in the partonic processes, 
while the hadronic jet with momentum Px contains all the hadrons as well as the beam remnants 
in the final state, except the top quark. Explicitly, px = Pi + k, where p\ is the momentum of the 
final state collinear partons forming the jet and k is the momentum of the soft radiations. Such 
division of momentum is artificial and we have to integrate over the soft momentum to obtain a 
physical observable. 

We also define the Mandelstam variables as 

s=(P a + P b ) 2 , u = (P a -q) 2 , t = (P b -q) 2 (19) 

for hadrons, and 

s = (p a +Pb) 2 , u = {Pa~q) 2 , i={p b -q) 2 (20) 
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for partons, respectively. In terms of the Mandelstam variables, the hadronic and partonic threshold 
variables are defined as 

S 4 = P x = s + t + u-m 2 , (21) 
s 4 = p 2 x = s + i + u — mf, (22) 



where mt is the mass of top quark. The hadronic threshold limit is defined as S4 — > [53]. In 
this limit, the final state radiations and beam remnants are highly suppressed, which leads to final 
states consisting of a top quark and a narrow jet, as well as the remaining soft radiations. Taking 
this limit requires x a — > 1, Xb — > 1, s 4 — > simultaneously. In this limit, we get 

S 4 = s 4 + »(— - 1) + (i - m 2 )(— - 1) + (u- m 2 t ){— - 1) 

X a Xb X\) x a 

« s 4 + s(x a + Xb) + m 2 )xb + {u — m 2 )x a 

~ «4 + (-i)x a + (~u)x b , (23) 

where x a ^ = 1 — x a ,b- This expression can help to check the factorization scale invariance which is 
shown with more detail in the following. The hadronic threshold enforces the partonic threshold. 
However, the reverse is not true. The partonic threshold s 4 — > does not forbid a significant 
amount of beam remnants. We note that in both hadronic and partonic threshold limits, the top 
quark is not forced to be produced at rest, i.e. it can have a large momentum. For later convenience, 
we can also write the threshold variable as 

S4 = Px = (Pa +Pb~ qf =p\ + 2k + E 1 + 0(k 2 ), (24) 

where k + = nrk, k is the sum of momenta of soft radiations; E\ is the energy of the quark jet and m 
is the lightlike vector associated with the jet direction. In the threshold limit (s 4 — > 0), incomplete 
cancelation between real and virtual corrections leads to singular distributions a™ [ln m (s 4 /m^ )/s 4 ]+, 
with m < 2n — 1. It is the purpose of threshold resummation to sum up these contributions to all 
orders in perturbation theory. 

The inclusive total cross section of the t-channel single top production can be written as 

1 da; 



a = J dx a J dx b j di j dufi/p^F^^fj/p^FjXb) 



'1.1 



2s dtdu 

I dp T / dy dx b / ds 4 — ■ 2\fi/P S^F^^fj/P b^'^^FJT > 

Jo J- ymax Jx b>min Jo 2(x b s + u-mf) dtdu 

(25) 
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where we have changed the integration variables to be the top quark transverse momentum squared 
Pt, rapidity y, Xb and S4. The regions of the integration variables are limited by 



PT,max 



4s 



1 \ + y/l-sq As(jp 2 T + m 2 t ) 

-in , with sq — 



2 H (s + m 2 ) 2 ' 



s + t — mf ' 



x fe (s + t - m£) + u, (26) 



i = m 2 — \fsJ p\ + m 2 e v 



with 



u = mf - VsyPr + m t e V - ( 27 ) 
The other kinematical variables can be expressed in terms of these four integration variables. 

IV. FACTORIZATION AND RESUMMATION FORMALISM 

To derive a factorization formula for the t-channel single top production in SCET, we first have 
to match the full QCD onto the effective theory. In this section, we follow the convention and 



formalism in [54, |55j, where the matching is performed in momentum space. The relevant operator 



in QCD responsible for the t-channel single top production is 

0(x) = (d^P L ui^P L b)(x), (28) 

where we have adopted the Feynman gauge for the W boson propagator. This operator contains 
three massless quarks, which can be described by collinear quarks in SCET, and a massive quark, 
which can be described by heavy quark effective theory |56fl . The presence of three different 
directions and a massive quark is a feature of single top production at hadron colliders. The 
operator 0(p), which is the Fourier transform of O(x), can be written in terms of the momentum- 
space SCET fields in the threshold region as 

n f d 4 p a d A p b d A px d 4 p 2 d 4 k s 
° {P) = J (2nf (2rr)< (2nf (27r) 4 C '(Pa,P 6 ;Pi,p 2 ) 

xO in (p a ,p b )0 ou t{p 1 ,p 2 )Os,i(k s )(2iT) A 6^{p - Pa - Pb +pi +p 2 + k s ), (29) 

where the operator Oi n (p a ,Pb) is responsible for annihilating the initial u- and b-quark with mo- 
menta p a and pb, respectively, or explicitly, 

C c$,in=XUPa)xi(p b ), (30) 
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and 0U t(pi,P2) is responsible for creating the final d- and t-quark with momenta p\ and p2, 
respectively, or explicitly, 



■y8,out 



Xs(Pl)KJP2 



(31) 



Note that we have described the top quark in terms of the heavy quark effective field with a 



label velocity v [56]. Since there are two fermion lines in this process, either of which connects an 
initial and an final states, we indicate them with the Lorentz (a, /3, 7, 5) and color indices (c, d, e, /) 
explicitly, and retain only quark fields in the operators for simplicity, leaving the other structure 
in the matching coefficient Ci, which is at the LO level 

g 2 V ud V tb 



(32) 



Here, the electroweak coupling is defined by g 2 = 8GfM^/^/2 where Gf is the Fermi constant. Vij 
is the CKM matrix element and Mm is the mass of the W boson. 6n denotes the color structure 
of the t-channel single top production at LO in the singlet-octet basis 



|ci)=5/ c 5ed, |C2> = (t a )fc(t a ) 



edi 



(33) 



where t a is the generator of the gauge group SU(3) C , satisfying 



tr[t a t b ] = -S ab . 



(34) 



I = 1 or 2 is an index in this color space. In Eq. (29), we can separate the soft gluon field from 
collinear or massive fields because of the field redefinition in Eq. (S) . 

The soft operators Os,i, which are responsible for the soft interactions between different collinear 
sectors and top quark, are expressed as 



O f s f{k s ) 



d 4 xe~ iks - x T 



d 4 xe~ iks - x T 



fc 



(Y,}{x)Y nb {x) 
Yj(x)t a Y nb (x) 



(35) 



where the time-ordering operator T is required to ensure the proper ordering of soft gluon fields 
in the soft Wilson line. 

The total cross section for t-channel single top production in the threshold region can be written 
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as: 



^(I\Ot(0)\X){X\O x (0)\I)(27r) 4 5HPa + P b -px) 

X 

±- s J2[ d*x(I\Ot(x)\X)(X\O x (0)\I) 



X 



= ^E/ ^j('|Ot(0)|X)(X|O(p)|/), (36) 

where |J) = \P a Pf,) denotes the initial state protons (anti-protons). Here we distinguish the position 
space operator from the momentum space one by a subscript x. The restriction on the sum over 
final states \X) is that the final state configuration consists only of a top quark jet whose 3- 
momentum is in the direction of n\, an d-quark jet in the direction of m, and soft radiations. This 
is the configuration that is relevant to threshold resummation and we are interested in. Under this 
condition the final state can be written as \X) = \XtX\X s ) , where \Xt), \X\) and \X S ) denote top 
quark jet, d-quark jet and remaining soft radiations, respectively. In the second line of Eq. (|36|), we 
have used the momentum conservation delta function to shift the operator O x to point x, and in 
the third line written the operators in momentum space, which are matched onto SCET operators. 
Using the notation <3?2 = {Pa,Pb]Pi,P2} to express a phase space point [54] with d$2 = 



d p a d Pbd pid p2 /(2tt) and $2 — k s = p a + Pb — Pi — P2 — k s , we can write Eq. (36) in a compact 
form 

- = ^E/ d^cK^c^) J |||^(2-)V 4 )(ci>2-^) 

x (/] (O'^O'^O's^ \XtX.Xs) {X t X x X s I (O in out O s ,i)\I) . (37) 

As we mentioned before, different collinear sectors are decoupled due to field redefinition, and thus 
the matrix element in Eq. ( |37| ) can be factorized into a product of several matrix elements, which 
obey certain RG equations. 

In the following, we further show the matrix elements mentioned above. First, we deal with the 
top quark sector. Since we have decoupled the soft interaction by field redefinition, the top quark 
now should be regarded as a non-interacting particle, which can be written as 



x t 
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where summation over final state \Xt) gives rise to a top quark phase space integral. Next, we 
define the soft function by the soft matrix element as 

(X s \O f s ^ d (k s )\0)5(k + - nvk a ), (39) 
where N c is the number of colors and we have inserted into the above equation an identity operator 

1 = J dk + 5[k + -m-ks], (40) 
because of the constraint from Eq. (24), which expresses the multipole expansion of a soft field 



interacting with a collinear field [H™. Note that the summation over final state can be performed 
Y2x a \Xs(k' s )) (X s (k s )\ = (27r) 4 5( 4 ) (k' s — k s ) since there is no restriction in the summation and also 
there is no explicit dependence of the final states on \X S ). Since we are only interested in the total 
cross sections, the final state top quark, jet function and PDFs can be considered to be diagonal 
in color space. Then we can contract their color indices to obtain the soft function matrix 

Sjj(k + ,fi) = 5f'f5 c ' c 5 e ' e 5 d ' d S d j? c ' f ' fced (k + ^). (41) 

At the LO, it can be written as 

j. 1 ( C\ \ 
S(k+,,) = 6(k+)-l * (42) 

where Ca is the Casimir operator for the adjoint representation of SU(3) C . At the NLO, the 
calculation of soft function boils down to the evaluation of eikonal diagrams p6| . Since the virtual 
corrections in SCET vanish, only real emission diagrams are needed to be evaluated. The details 
of the calculation of these diagrams are given in Appendix 
For the final state d-quark jet sector, we have 

j4„/ j4„ 



where the summation over collinear state has been performed and J is the spin- and color-singlet 
jet function, defined as 



d(p°)n 1 .pJ(p 2 ) = ^ W J j^Tr(0\rjt lX (p')x(p)\0), (44) 
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where TV represents the trace over spin and color indices. At LO, it is just 5(p 2 ). Finally, the 
initial state n a collinear sector reduces to the conventional PDFs: 

d A p' a d 4 p a 



(Pa\x C a'(Pa)Xa(Pa)\Pa) 



JL 5 c'c f 1 ^ L aEc JA f u/Pa (x a ,v) (45) 



(2vr) 4 (2^" awya ^ a ^ aKian a ' 2N C 

and similarly for the matrix element for n b direction. Thus the momenta of incoming partons are 
given by p a(b) = x a(b) EcMn a {b)/2. 

Combining the above expressions, we obtain (up to power corrections) 

„ 1 daj^ 
dx a dx b dtdu—fi/p a (x a , n)f j/P {x b , y) / , (46) 
2s ' dtdu 



a 



with 



da^ cs 1 1 1 



and 



dtdu 4iV c 2 8vr s 

x f dk + J dp 2 1 Sj I (k + , f i)J(piv)5(s 4 -p 2 1 -2k + E 1 ), (47) 



\ov = g 2 \v ld \ 2 \Vjt\ 2 S m pt - (48) 



All the objects in the factorized Eq. ([46]) have precise field-theoretic definitions so that they can be 
calculated directly and systematically, except the nonperturbative PDF. The convolution between 
the jet and soft function suggests that the partonic threshold consists of two parts. In the case 
of S4 = 0, there are no collinear or soft gluons emitted. In the small S4 region, the number and 
momentum of collinear and soft gluons are constrained. 

At the LO, the hard function Hjj is normalized to Sh5ji. In general, it is related to the 
amplitudes of full theory, which is given by [^] 

X 04j HQ = —-L-— ({ Cl \M<&>(Mj®\cj) + (ci\M£l}(M£l\cj}) , (49) 
\ci\ci){cj\cj) V / 

where (A'lren) are obtained by subtracting the IR divergences in the MS scheme from the UV 
renormalized amplitudes of full theory. 

Because of the special color structure of this process, the hard function matrix elements do not 
contribute to the cross section except for H\\ at the NLO level. In SCET, there is a RG evolution 
factor connecting the hard scale [ih and the final common scale //, which would contain contributions 
from non-diagonal elements beyond NLO. However, these non-diagonal contributions involve the 
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gluon connecting two fermion lines, resulting in a suppressed color factor 1/N 2 , compared to 
diagonal ones. Thus, we expect their contributions are small and can be neglected safely. Then 
the t-channel single top production is considered to be a double deep- inelastic-scattering (DDIS) 
process ||. In this case the hard function H%\ can be further factorized into two parts, i.e., H up 
and Hdn, which represent contributions from the up and down fermion line in Feynman diagrams 
as shown in Fig. [I], respectively. This separation is also helpful to make a reliable perturbative 
prediction for the hard function. The reason is that usually the loop corrections from the up 
and down fermion lines contain large logarithms of the form ln(— t/fj,V) and ln((— t + mf)/mi//%), 
respectively; see Eqs. (|5"l|)-(|5"2|). It is hard to choose a proper hard scale to make both of them 
small. In the case of a DDIS process, the two separate hard parts can be evaluated in different 
scales such that the perturbative expansion is reliable in both parts. As a consequence, we can 
rewrite the Eq. fl4"7|) as 

dajf res l 1 1 

x J dk + J dp 2 l S(k + ,fi)J(p 2 1 ,fj 1 )5(s 4 -p 2 l -2k + E 1 ), (50) 
where S(k + ,ii) denotes the component Sn(k + , fi) in Eq. ([4l|). 

V. THE HARD, JET AND SOFT FUNCTIONS AT NLO 

The hard, jet and soft functions describe interactions at different scales, and they can be calcu- 
lated order by order in perturbative theory at each scale. At the next-to-next-to-leading logarithmic 
accuracy, we need the explicit expressions of hard, jet and soft functions up to NLO. 



A. Hard functions 



The hard functions are absolute value squared of the Wilson coefficients of the operators, which 
can be obtained by matching the full theory onto SCET. In practice, we need to calculate the one- 
loop on-shell Feynman diagrams of this process in both the full theory and SCET. In dimensional 
regularization (DR), the facts that the IR structure of the full theory and the effective theory is 
identical and the on-shell integrals are scaleless and vanish in SCET imply that the IR divergence 
of the full theory is just the negative of the UV divergence of SCET. After calculating the one-loop 
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(a) 




FIG. 2: The one-loop Feynman diagrams for the t-channel single top production. 

on-shell Feynman diagrams, as shown in Fig. we get the hard functions at NLO as follows, 

H up (» hiUp ) = 1 + Cpa j^ f-21n 2 ^ + 61n^- + , (51) 

47r V Vh,up Vh,up ) 

tt i \ C F a s {n h4n ) ( 2 -i + m 2 t -i + m 2 t H dn \ 

Hdn{fJ-h,dn) = H : ! — -41n z h lOln h c x , (52) 

4vr V Vh,dnm t Hh,d n m t J 



where 



2 



cf'" P = -16 + y, (53) 

cf > dn = ~ ln(l - A) + 21n 2 (l - A) + 6 ln(l - A) + 4Li 2 (A) - 12 - — 
A 6 

2m 2 u . m 2 . A . 

+ tr~ 2> ^ ± r ( 54 ) 
t{s — mi) mf — t 

with A = t/(i — m 2 ). These results agree with those in Ref. Q. In order to avoid large logarithms, 
the natural choices of Hh,up and ^h,dn are \f— i and (—i+m 2 )/m t , respectively. The hard functions 
at the other scales can be obtained by evolution of RG equations. The RG equations for hard func- 
tions are governed by the anomalous-dimension matrix, the structure of which has been predicted 



up to two- loop and four- loop level for the case involving massive [57] and massless partons 
respectively. In our case, we can write the RG equations for hard functions as 

j-j H up (n h ^ up ) = ( 2r cusp In— h 2jV p J H up (fih,u P )i (55) 

r^h,up V ^h,up ' 

d f —t + m 2 



£ _|_ jyi^ \ 

2r cusp In + 27^ H dn (n h dn ), (56) 

^h,drjn t J 



where r cusp is related to the cusp anomalous dimension of Wilson loops with light-like segments [59], 
while 7^1 and j dn control the single-logarithmic evolution. Their expressions up to two-loop level 
are shown in Appendix [B]. 

After solving the RG equations, we get the hard functions at an arbitrary scale /i: 

H U p(/J,) = exp[AS(/j,h, U p,IJ) ~ 2 a«p(M/i,«p>M)] ( — ) H up (/j, hiUp ), (57) 

Hdn{y) = exp[25(/i Mn ,/i) - 2a^ n (fi hjdn , fi)] ( H dn (/j, h dn ), (58) 
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where S(nh tUp ,fi) and denned as |6(| 



S(/J>h,up, M) 



"'up 



a s (» -p / 

, 1 cusp V 



a) /*" da' 

) A.(A. h ,«,) /%') 



a s (tJ, h>up ) A%) 



(59) 

(60) 



S{jJ>h,dni A 4 ) j a r and have similar expressions. 



B. Jet function 



The jet function J(p 2 ,/u), defined in Eq. (f43|), describes a quark jet of invariant mass squared 
p 2 . It is process-independent and has been calculated at NLO in and NNLO in 61]. The RG 
evolution of the jet function is given by 

^f^ = (- 2 r ra!p i„4- 2 y) J(A«) + 2r cmp r'tf ^-rt-y.rt . (61) 

a In// V A* ) Jo P ~ Q 

This integro-differential evolution equation can be solved by using the Laplace transformed jet 



function [25, 6C]: 

12 r°° J2 



~ Q 2 f 00 p 2 

j(ln— ,//)=/ dp 2 exp(- )J(p 2 ,fi), (62) 

A* Jo y z eTs 



which satisfies the the RG equation 



j(ln%^)= (-2T cusp ln%-2 1 J )j(ln%, f ,). (63) 



a in // // z \ /i z J fj, z 

Then the jet function at an arbitrary scale // is given by 

~ 1 / p 2 \ e~ lEll: > 

J(p ,//) =exp[-4J5(jij,n) + 2a (mj>aO]j(\>Mj)^ ( J£ J "r^p ^ 

where r/j = 2ar(//j,A t )- The /i-dependent part of the Laplace transformed jet function j(L,//) is 
determined by the anomalous-dimensions of the jet function as in Eq. (^), while the //- independent 
part can only be obtained by a fixed-order calculation. At NLO, it is 

a x /Tn 



j{L,n) = l + ^-l^L 2 + ^L + ci), (65) 



with c{ = (7- §7T 2 ) Cf. 
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FIG. 3: Nonvanishing diagrams contributing to the soft function at NLO. The contributions from the left 



and right diagrams are denoted as Sju and , respectively. 



C. Soft function 



The soft function S(k + ,Li), denned in Eq. (|39"|), describe soft interactions between all colored 
particles. It can be calculated in SCET or in the full theory in the Eikonal approximation. The 
LO soft function is given in Eq. ([i2|). At NLO, we only need to calculate the nonvanishing real 
emission diagrams in dimensional regularization, as shown in Fig. ||, which give 
(1) 2gjC F ^ 1 [°° + f°° _f , 2 

and 



q±S(q + q - q±)6(k - n x ■ q) 



rib ■ v 



(q ■ n b )(q ■ v) 



(66) 



1 



(67) 



'0 JO J 

respectively. After calculating these integrals by the approach of Ref. 



we get 



47T 



In- 



+ 5{k)c s bt 



and 



S tt (k,n) = 5(k) 



2C F a s 
Air 





"2" 




{- 


k_ 


+ 5(k)4j 



(68) 



respectively, where jl = fi/y (2n h i) / 'nf 2 = (Lt(—u)m t )/(2(—i + mf)Ei). The detail of our calcu- 
lations and explicit expressions of cf t and cf t are given in Appendix B. The star distribution is 
defined as [39] 

[f(x)}^ = f(x) for x > 



dx[f{x)t^g{x) 



dxf{x) [g{x) - g(0)] . 



(69) 
(70) 



/o Jo 

And the soft function S(k,fi) = Sbt(k,Li) + Su(k,Li), similar to the jet function, satisfies the RG 

equation 

d 



din \i 



S(k,Li) 



2r cusp ln- + 2 7 - 



S(k,fi)+2T. 



cusp 



k ^ S(k^)-S{k'^) 

dk J—kl • 



(71) 
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The solution to this equation is 



S(k,/j) = exp[-25(/i a ,/i) - 2a s (fi s ,fi)]s(d Va ,fi s )- 



1 / Jfc \ Va e~^ s 



k\£i s J T(r] s ) ' 

where rj s = 2ar(/u s ,/u). The Laplace transformed soft function s(L,/i) at NLO is given by 



with cf = (2cg + 2c| + ^-)C F 



(72) 



(73) 



D. Scale independence 



In the factorization formalism, we have introduced the hard function, jet function and soft 
function. Each of them is evaluated at a scale to make the perturbatvie expansion reliable, and 
then evolved to a common scale. Therefore, it is important to check the scale independence of 
the final results. If we expand the exponent in Eq. (|57|), then we can find the dependence on the 
intermediate scale fJ,h,up cancel each other up to 0(a s ). The same situation happens for Hh,dn in 
Eq. (|5^). The case for the jet scale is more complicate due to the partial derivative operator and 
the delta function after we use the expansion 



1 



•h 



p 2 u? 



6(p 2 ) 



+ 



1 

p 2 



+ Vj 



ln(p 2 /n 



P 



(74) 



We point out that the scale independence happens for the jet function only in the sense of the 
integration over p 2 . The case for the soft scale is the same as for the jet scale. 

After checking the intermediate scale independence, we discuss the case for the final common 



scale. Recalling the hadronic threshold definition in Eq. (23) and the cross section near the 
threshold in Eq. (|50|), we have 
da 



dS^dy 



dx a dx b J dpj J dk + jf i/Pa (x a , (i)f j/Pb (x b , n)H up {fj)H dn {n) 
J(p 2 , f i)S(k+, f i)5(S i - (-£)(! -x a ) - (-u)(l-x b )-p 2 -2k + E l ), 



(75) 



where we have changed the integration variables dtdu to dp^dy and then to dS^dy. From this 
equation, we can see clearly the connection between the threshold region of the whole system, 
represented by S4, and those of the parts of the system, represented by (1 — x a ), (1 — x b ),p\ and 
k + , respectively. To change the convolution form to simpler product form, we apply the Laplace 
transformation to the above equation and obtain 

do f°° ,„ f S A \ da 



dQ 2 dy 



dS - exp , . I • 

o V Q en E J dS A dy 



(76) 
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The Laplace transformed jet function and its RG evolution are given in Eq. (^) and Eq. 
Here, for convenience, we write its RG equation again as 



d ~ Q 2 

— J (In— 

rfln/j ^ 



Q 2 



-2r„, sn In ^ - 2 7 J ) j(ln ^, fi). 



CUSp n 



Q 2 



(77) 



The Laplace transformed soft function is similar to the jet function, but its RG equation is 



d Q^-j+m 2 ) 
-s(ln — ,n) 



dln/i v n(—u)mt 
The Laplace transformed PDF near the end point is given by 

1 - x 



J n(-u)mt 



2r cuS p In 

H{-u)mt 



fi/p(r, At) 



dx exp 



fi/p a (x,v), 



A*)- (78) 



(79) 



which satisfies RG equation 

d 



dlnfi' 



f i/P (T,fj) = 2r cusp ln(r) + 2^\ f i/P (r,n). 



(80) 



Due to the delta function in Eq. (|75|), the variables r a & in the Laplace transformed PDF are given 
by 



for 



fi/p a ( T a,v), and r b 



Q' 2 



for f j/Pb (n,^. 



(81) 



Q 2 

T a = —z 

—t —u 
For completeness, we also need the RG equations for the hard functions which have been given by 
Eq. (|55| ) and Eq. ([[(j). We rewritten them as 



d 



d ln/i 
d 



H up (n) 

Hdnip) 



-t 



2r cusp In + 27^, JflwpC/i), 



2r cuS p ln- 



-i + mf 



+ 2^ n )HM. 



(82) 
(83) 



d ln/i \ //mt 

So far, we can check the scale independence of the final results. Using the relation between 
anomalous dimensions given in Eq. (|B6| ), we can immediately obtain 
d 



dlnfi 

Even more precisely, we have 



fi/p a ( T a, v)fj/p b (n, fj)H up (/j,)H dn (iJ,)j(]xL^,n)s(]n - \ t ^ mt \ ^i) 
1 Jl fi z fjL{—u)mt 



0. 



(84) 



din fi 



~ Q 2 
fi/P a ( T a, v)Hu P (v)j(ln — ,fi) 

I- 1 



d 



dlnfi 



Q 2 



-t + mi 



>M) 



0, 
0. 



(85) 



H(-u)m t 

It means that if we evolve the scales of H up and the jet function to the factorization scale of the 
light quark line fiF,up, then final results should not depend on fiF,up- And if we evolve the scales of 
Hdn and the soft function to the factorization scale of the heavy quark line fiF,dn, then final results 
should not depend on fiF,dn- Actually, the relationships between the anomalous dimensions given 
in Eq. (|B6|) are determined by these requirements. 
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E. Final RG improved differential cross section 

After combining the hard, jet and soft function together, we obtain the resummed differential 
cross section for t-channel single top production 



dtdu ^ 64vriV c 2 s 2 



exp[4S(fj> h>up , hf,u P ) - 2a up (fi h:Up , fiF,u P )] ( — — ] H up (fj, htUp ) 

(—1 _|_ m 2 \ - 2a r(Hh,d n ,VF,d n ) 
Hdn{^h,dn) 
l^h,dnm t J 



exp[-4S(fij, fi F , up ) + 2a J {fij, fi F , up )] ( ^rj 



exp[-2S(fj, s ,fj, Fidn ) - 2a s {n s ,fi F , dn )] ^ 



m t {—i + mf t 
/i s (-it) 



'Is 



j(d v + L j ,fi j )s(d v + L s ,fi s )—[ — \ — , (86) 

where r] = r]j+r) s and Lj = ln(m 2 //i 2 ), L s = ln[mt(— t + m 2 )//x s (— u)]. In the above expression, the 
hard function H up and jet function {H dn and soft function) have been evolved to the scale fiF,up 
{l^F,dn}- It seems that the t-channel single top production is factorized as two DIS processes. How- 
ever, the convolution of the jet and soft functions, now expressed in terms of the partial derivative 
operator acting on the same kernel function, violates this simple factorization and connects the 
two DIS processes nontrivially. 

If we set scales Hh,up, V-h,dm t 1 -]-, Ms equal to the common scale //, which is conveniently chosen as 
the factorization scale fi F;U p = V>F,dn = fJ-F, then we recover the threshold singular plus distributions, 
which should appear in the fixed-order calculation. Up to order a 2 , we have 



+ [ B ^ + B ^ + B ^ + B i D i + B oH^, (87) 



where 



D n = 



ln n - 1 (s 4 M 2 ) 
s 4 



(88) 
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The A n and B n coefficients are given by 

A 2 = 3r , (89) 
Ax = (^ + 2^0+70-270, (90) 

Ao = ^-^Mp ~~ L h,dn + + L l- r - lup,0 L h,up ~ 2jdn,0 L h,dn + 76%' ~ 2jo L s 

+cf + cf + cf , (91) 
= (92) 

5 3 = ^Air - ^ r , (93) 

37T 2 

5 2 = A 2 + A 2 A - — - /3 r (L i + 4L S ) - /3 (7 J - 4 7o s ) + 3T l9 (94) 

Bi = Ai (a - yr ) + 9C 3 rg - /3 r QlJ + 2L 2 - Ayr 2 ) - /3 (cf + 2cf + 7 %- - 4 7o %) 
+7i J -2 7 f + T l (L 3 +2L s ), (95) 

where Cs = 1.20206 • • • , L h>dn = ln[(-t + m 2 )/ j Ltm t ], L Mp = ln(-t7/x 2 ) and cf = cf >ttp + cf ' dn . We 



find that A2 i i, -64,3,2 and the scale-dependent parts of Aq and £>i agree with the results in Ref. [2C] 
with the replacement t(u) — > t(u) — m 2 due to the different definition of S4 there. 

To give precise predictions, we resum the singular terms to all orders and include the non- 
singular terms up to NLO. We obtain the final RG improved differential cross section 

d&Resum _ dfi tlalCB ( d&^Q _ dtf* 



(96) 

expanded to NLO 



dtdii dtdu V dtdu dtdu 

Near the threshold regions, the expansion of the resummed result approaches the fixed-order one 
so that the second term in the above equation almost vanishes and the threshold contribution 
dominates. In the regions far from the threshold limit, the resummation effect is not important 
and the final result is mainly determined by the fixed-order calculations. 

VI. NUMERICAL DISCUSSION 

In this section, we discuss the numerical results for threshold resummation effects on t-channel 
single top production at the Tevatron (-y/i = 1.96 TeV) and the LHC (yfs = 8 TeV). The top quark 



mass is chosen as mt = 173.2 GeV [62] and the rapidity is integrated over —2.4 < y < 2.4 if not 



specified explicitly. For the W boson mass we take My/ = 80.4 GeV. We set the Fermi constant 
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FIG. 4: The Feynman diagrams for the single top production via the fusion of a W boson and a bottom 
quark. 

to be Gf = 1.1664 x 1(T 5 GeV" 2 . The CKM matrix is given by 

( 0.9751 0.2215 0.0035 ^ 
Vckm = 0.2210 0.9743 0.0410 (97) 

V 1 J 

Throughout the numerical calculations, we use the MSTW2008nnlo PDF sets and associated strong 
coupling constant. The factorization scales are set at m< unless specified otherwise. There are four 
other scales, i.e., fih,up, Hh,dn, fJ>j, [J-s, introduced in the factorization formalism. They should be 
properly chosen so that the corresponding hard functions, jet function and soft function have 
stable perturbative expansions. In order to achieve this aim, each function should better not 
contain large logarithms. From Eqs. (|5ll)-(|52]), we can see that if we choose Hhuv = Q = v~i an d 
l*>h,dn = {Q 2 + fn 2 )/^, then the large logarithms disappear. Also as discussed below Eq. ([49|), 
if we combine the two hard functions blindly, we can not choose a proper hard scale to eliminate 
all the large logarithms simultaneously. This is due to the fact that intrinsically the W boson 
connects interactions at different scales. We can take another viewpoint on the t-channel single 
top production and consider it to be a fusion process, as shown in Fig. An initial-state up 
quark emits a W boson, which then combines with a bottom quark to produce a single top quark. 
The production of the W boson is similar to a DIS process and there is no specific constrain on 
the virtuality of the W boson. But when it co-annihilates with a bottom quark, the mass of final- 
state top quark impose constraints on the 'initial' W boson. As a result, the typical scales of the 
interactions involving the light quarks and top quarks are Q and mt, respectively, which are just 
about the natural hard scales. 

For the jet and soft scales, the situations are not so clear. After inspection of Eqs. (|6^)-(|65|) and 
(|72])-([73D, one finds that the natural jet and soft scales should be \fp* and 2kEi(—i+mf)/(—u)/m t , 
respectively. But these two scales are not directly connect to the integration variables in Eq. (| 
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FIG. 5: The contributions from jet and soft functions to the NLO cross section. 
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W.up/Q W,,dn/((fi 2 +m» 2 )M) 
FIG. 6: The contributions from hard functions to the NLO cross section. 

Moreover, they can become so small that strong coupling constant in the jet and soft functions 
would diverges. Therefore, in practice, we choose the natural jet and soft scales numerically. In 
Fig. ||, we show the contributions to the NLO cross section from jet and soft functions separately 
without including the RG evolution effects. We have fixed the top quark transverse momentum 
to be pr = 20, 60, 100 GeV and change the jet (soft) scale from 5 GeV to 100 GeV. It is required 
that the perturbative expansions of the jet and soft functions converge fast. Thus, we choose the 
jet and soft scale as 80 GeV and 50 GeV, respectively. When giving the final RG improved cross 
sections, we will investigate the scale uncertainties due to these choices. From Fig. ||, we can also 
see that the jet and soft function give positive contributions to the NLO cross sections, and can 
be as large as about 20% and 40%. To see the corrections from hard functions, in Fig. ^, we show 
the contributions to the NLO cross section from hard functions. We find that the hard functions 
provide negative contributions to the NLO cross sections and the corrections are about —15% for 
Vh,u P = Q and fj, hjdn = (Q 2 + m 2 )/m t . 

Before presenting the numerical results for RG improved cross section, it is important to examine 
to what extent the singular terms approximate the fixed-order calculation. In Fig. we present 
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FIG. 7: The singular terms contribution and fixed-order cross sections for t-channel single top production 
at the Tevatron (left) and the LHC (right). The dashed line represents the contributions from the singular 
terms up to 0(a s ) which is given in Eq. (p7|). 



the singular terms contribution and fixed-order cross sections. We see that the NLO cross section 
is well approximated by the singular terms when the top quark transverse momentum pt is larger 
than 50 (70) GeV at the Tevatron (LHC). Therefore, the singular terms should be resummed for the 
large pt region. For the small pr region, the singular terms do not dominate the NLO corrections, 
so there is no need to perform resummation in this region. In the following discussion, we will only 
present the resummation results for pr > 50 (70) GeV at the Tevatron (LHC). Meanwhile we find 
that the NLO QCD correction is small for t-channel single top production. This is because the 
large positive soft and jet functions cancel with the large negative hard functions, as discussed in 
the last paragraph. If these large effects are resummed to higher orders, we can see whether there 
is still a cancelation between them. 

We have chosen all the natural scales involved in this process. Now we give the numerical 
results of the resummed cross section. When discussing each scale dependence, we fix the other 
scales at the natural scales discussed above. In Fig. H we show the RG improved cross sections as 
a function of the top quark pr- We can see that the distribution is increased by about 9 ~ 13% 
and 4 ~ 9% for pj- > 50 and 70 GeV at the Tevatron and LHC, respectively, compared to the 
NLO results. In Fig. §, we give the uncertainties of the resummation results due to the change of 
intermediate scales fJ<h,up-> ^h,dni Mj> Ms independently by a factor of two. The uncertainties arising 
from Hh,up, ^h,dn and fij are less than ±1%, and for fi s are about ±2%. In Fig. 1C, we show the 



scale uncertainties of the resummation results due to the variations of HF,dn and /iF,up by a factor 
of two, and do not see scale uncertainties are decreased, compared to the NLO results. In principle, 
the scale uncertainties should vanish, as illustrated analytically in the last section. However, the 
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FIG. 8: The RG improved (dashed) and fixed-order cross sections for t-channel single top production at the 
Tevatron (left) and the LHC (right). We have defined <5crR CSUm = dan csum — dcNLO- 




FIG. 9: The scale uncertainties of the resummation results due to the variations of Hh,upi ^h,dn, V>j, Ms, 
respectively. 
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FIG. 10: The scale uncertainties due to the variations of (j*f dn (left) and PF,up (right), respectively. 

analysis there is based on the assumption that the PDF is evaluated near the endpoint. But in 
practice, this is not always true because the center-of-mass energy of the Tevatron or LHC is much 
larger than the invariant mass of the final states. And the dynamical enhancement mechanism [28] 
is not appropriate for a t-channel process. On the other hand, when approaching the threshold 
region, i.e., with the increasing of the top quark pr, the scale uncertainties of the resummed cross 



sections are significantly reduced, as shown in Fig. 10 



VII. CONCLUSION 



We have studied the factorization and resummation of t-channel top quark transverse momen- 
tum distribution at large px m the SM at both the Tevatron and the LHC with SCET. This is the 
first space-like process studied in SCET involving one massless and one massive colored particles 
in the final states. The cross section in the threshold region can be factorized into a convolution of 
hard, jet and soft functions. In particular, we first calculate the NLO soft functions for this process, 
and give a RG improved cross section by evolving the different functions to a common scale. Our 
results show that the resummation effects increase the NLO results by about 9 ~ 13% and 4 ~ 9% 
when the transverse momentum of the top quark is larger than 50 and 70 GeV at the Tevatron and 
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the 8 TeV LHC, respectively. Our prediction on the transverse momentum distribution of the top 
quark in the large pt region is important in the searching for new physics, e.g., a heavy W' which 
can mediate the single top production through s-channel. Also, we discuss the scale independence 
of the cross section analytically and show how to choose the proper scales at which the perturbative 
expansion can converge fast. 
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Appendix A: Calculation of the soft functions 

In this Appendix, we present the details of the calculation of the two 0(a s ) soft functions 
5^(/c,/i) and (k,fJL). We choose to do the calculation in the rest frame of top quark, in 
which the four-velocity of the top quark is = (1,0,0,0). This choice of the frame makes the 
denominators simple but leave the complexity in the delta functions. Actually, we also perform the 
calculation in the frame where the delta functions is simple but the singularities in the denominators 
are hard to isolate [53]. And finally we find the same results, which can be considered as a strong 



cross check for our calculations. 

In rest frame of top quark, we also choose = (1,0, 0, 1). Then 



and 



q +< +q -< +q ^ n V = n +< +n -< +n ^ (A1) 
n bb n bb n bb n bb 



q + n 1 +q nf . ., . (n b + n- b ) (q + + q ) 
q-ni = — cost', q-v = q- — — = . (A2) 



■n 



hi, 



After putting these expressions into the formula (pq), we get 



xn q n i +q n i i i n i a\ nb ' v / a o\ 

5{k ~ + \q±\\ni±\cos6)——— -. (A3) 

n hl q + {q + + q-) 



Now redefine the integration variables q + and q and let a = -h, then 

n i 

J2r~f ,,2e roo poo 



6{k -q+- q -+ 2^qT q ^ COS e)-^^— ) . (A4) 
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Introducing two variables x and y such that q + = kyx and q = ky(l — x) = kyx 
Jin ,,2e 



(27T) 



jfe" 



The singularity in the integrand can be isolated by 



-1-e 



1 







lnx 



ax + x 



After completing the above three parts of the integration separately and expanding 



2c 



2e 



<5(fc + ) + 



[*+,£] 



2e 



1 fc" 
y-ln- 



[k+,p\ 



+ 0(e 2 ), 



we get the divergent and finite parts 



2C F a 8 {4nrn 2 e-^ E Y ( 6{k) 2 



4tt 



V 



n Ml 



2C F a s 
4tt 





rln^i 











+ <5(A:)c; 



withcf i = -ln 2 (l + I)-2Li 2 ( T ^) + fJ. 
In the same method, we can get 



4^ {"Tj' 



2C F a s 

47T 





to 


{- 


A 



[fe,/i] 



(A5) 



(A6) 



(A7) 



(A8) 
(A9) 



(A10) 
(All) 



with 4 = 21n(l + A). 

When performing the Laplace transformation from S(k,fi) to s(L,fj,), we use the following 
replacements: 



In * 1 IM 



[M 



L 2 vr 2 
2 12' 



L. 



(A12) 
(A13) 



Appendix B: anomalous dimensions 



The various anomalous dimensions needed in our calculations can be found, e.g., in [28, 46, 60]. 
We list them below for convenience of the reader. The QCD /3 function is 

2 



P{<Xs 



-2a, 



(Bl) 
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with expansion coefficients 



ft 



2857 



Ci + ( 2C 2 F 



2 ^c F c A 1415 ^ 



44 



c l ) T ? n f + ( if + iff <^ ) I$nJ, 



158 
~27~ 



54 " V 9 27 

where = 3, C F = 4/3, Tp = 1/2 for QCD, and n/ is the number of active quark flavor. 
The cusp anomalous dimension is 

r cusp (a s ) = r ^ + r 1 (J) +••• , 



(B2) 



(B3) 



with 



r 
r. 



AC F . 
AC F 

AC F 



245 134 
~6 27~ 

55 



2 11 4 22 \ ^ ^ / 418 40 2 56 \ 

vr 2 H 7r 4 H (3 + C A T F n f 1 ir 2 C3 

45 3 7 * V 27 27 3 7 



+C F T F n f [ ~ + I6C3 



27 F y 



(B4) 



The other anomalous dimensions are expanded as Eq. (B3), and their expansion coefficients are 



7 ? ° 


= — 3(7f, 




= cl{- 




= -2C F , 


7Q 


= C f Ca 


7^ 


= 3Cp, 


7^ 


-<*(f 


7? 


= — 3Cp, 


7) 


= c|(- 



961 11 2 



- + 27T 2 - 24C 3 ) + C F C A ( -— - + 26C 3 ) + C^n, ( — + ) , 



4Cs ) + Y CFTpnf > 



2tt 2 + 24C 3 + C F C A 



17 22 



+ 2vr 2 - 24C 3 + C F C A 



242 



6 9 

1769 11 



54 



vr 2 + 40C; 



+C F T F n f [ — + -if 



130 



27 



+ -12C3)-C F T F n / (- + -7r' 



(B5) 



7up> 7dn ano - 7^ can ^ e obtained from the anomalous dimensions above through the following 
equations: 



1 V 
1 up 



27 g , 

7dn = Iq + lQ 
s 



V 



7 



70 -Idn 



(B6) 
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